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We study charge fluctuations within the statistical hadronization model. Considering both the 
particle yield ratios and the charge fluctuations we show that it is possible to differentiate between 
chemical equilibrium and non-equilibrium freeze-out conditions. As an example of the procedure 
we show quantitatively how the relative yield ratio A/A” together with the normalized net charge 
fluctuation v{Q) — (AQ^) / (Wh) constrain the chemical conditions at freeze-out. We also discuss 
the influence of the limited detector acceptance on fluctuation measurements, and show how this 
can be accounted for within a quantitative analysis. 
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I. INTRODUCTION 

In relativistic heavy ion collisions a localized high en¬ 
ergy density domain, a fireball, is created. The study 
of the properties of this hot and dense matter is the 
main objective of the experiments being conducted at 
RHIC and as of 2007 at LHC. Event-by-event particle 
fluctuations are the observables sub_^ct to intense cur¬ 
rent theoretical 0, S S S IE IE 0 IE IE i experimen¬ 

tal [IE in EE interest. Fluctuation measurements are 
important since they can be used: (i) as a consistency 
check for existing models, e.g. within statistical parti¬ 
cle production models 0 IE : (h) as a way to search for 
new physics, including C^P 0 m EE (iii) as a test of 
particle equilibration 00- 

The statistical hadronization model (SHM), intro¬ 
duced by Fermi in 1950 EEEEE3 , has been used exten¬ 
sively in recent years in the study of strongly interacting 
particle production. In this model, the properties of the 
final state particles are determined by requiring that the 
final state maximizes entropy given the physical proper¬ 
ties of the fireball (energy, baryon content, etc.). When 
the full spectrum of hadronic resonances is included EE, 
the SHM turns into a quantitative model capable of de¬ 
scribing in detail the abundances of all hadronic particles. 

Fluctuations in conserved quantum numbers (such as 
charge, baryon number, strangeness, or equivalently the 
net multiplicities of up, down and strange quarks) can 
be studied only in the Grand Canonical (GC) ensem¬ 
ble, since in the micro-canonical and canonical ensembles 
these quantities are fixed. We also mention here that fluc¬ 
tuations of non-conserved observables, e.g. other hadron 
multiplicities, differ for different ensembles even in the 
thermodynamic limit ^EiEE- 

In this paper we will discuss the use of fluctuations 
as a phenomenological tool within the framework of the 


statistical model, and illustrate some issues pertinent in 
analyzing fluctuations data. In section EH we will mo¬ 
tivate the choice of charge fluctuations as a useful ex¬ 
perimental probe. After demonstrating, in section EH 
how the statistical model implies a scaling between fluc¬ 
tuations and yields, we show (section E3 that a mea¬ 
surement of both particle yields and charge fluctuations 
can distinguish between an equilibrium high temperature 
statistical freeze-out from a super-cooled over-saturated 
freeze-out from a high entropy phase. Finally, in sec¬ 
tion El we discuss issues related to detector acceptance 
which impact the fluctuation measurement even in a 
boost-invariant azimuthally symmetric limit. We quanti¬ 
tatively demonstrate how such limited acceptance effects 
can be taken into account and the freeze-out temperature 
and non-equilibrium parameters extracted from experi¬ 
mental data. 


II. GC OBSERVABLES 

A study of GC SHM fluctuations of conserved quanti¬ 
ties is of considerable interest at RHIC. Since the detec¬ 
tors at RHIC (except for the PHOBOS detector) only see 
small portions of the final phase space, using the grand- 
canonical approach is justified in the following sense: 
Provided the fireball is indeed locally thermalized, we 
can take the experimentally observed source to be a sub¬ 
system in contact with a larger reservoir. 

The situation is of particular interest for reac¬ 
tions at RHIC that exhibit a sizable central plateau 
in the (pseudo-)rapidity spectrum, since a limited 
(pseudo)rapidity acceptance window selects a suitable 
subset of the source particles. Specifically, it can be 
shown (sections |^. The reasoning used there can be 
generalized to Fermi-Dirac and Bose-Einstein statistics) 
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that the rapidity spectrum of a boost invariant system 
could be related to the multiplicity in a static GC sys¬ 
tem with the same temperature and chemical potentials 

where i and j are species labels and the subscripts b.i. 
and GC denote the boost invariant system and the grand 
canonical system, respectively. 

The derivation in |2ll| can be applied to fluctuations at 
hadronization (before resonance decays) to show 


where Qi is the particle charge and i labels all particles 
before resonance decays since net charge is conserved Q . 

To use Eq. © quantitatively, however, the experimen¬ 
tal rapidity window must be large enough to encompass 
all decay particles of the resonances, yet small enough for 
the GC ensemble to maintain it’s validity. See section Fvl 
for a discussion of the validity of this assumption, and 
how to incorporate deviations from it in realistic experi¬ 
mental situations. 


III. STATISTICAL HADRONIZATION 


{dANi/dy\^ ^ 

{dN,/dy)^,^ 

where we denote the variance (fluctuation) of any quan¬ 
tity A as (AA2) = (a2) - {Xf. 

Given this, SHM average yields and yield fluctuations 
can be calculated by a textbook method , as per sec¬ 
tion ^1 

When studying finite systems the consideration of fluc¬ 
tuations in extensive quantities such as of particle yield 
has to address also volume fluctuations when the volume 
cannot be flxed by experimental conditions. In our case 
volume fluctuations can arise due to initial reaction ef¬ 
fects, impact parameter variations, as well as from fluc¬ 
tuations due to dynamics of the expanding fireball. It 
is difficult to arrive at a reliable description of all these 
effects. Therefore it is important to select fluctuation 
observables in which volume fluctuation effects are sub¬ 
dominant. Among extensive quantities, the net charge 
fluctuation stands out as it is relatively easy to measure 
and can be shown to be nearly independent of the volume 
fluctuations P|. 

In light of the above considerations we concentrate our 
effort on the following net charge fluctuation measure: 


For a hadron with an energy Ep = \/p^ + 'm?, the GG 
partition function for each species is given by 

In Z, = (t)E5. I -0^ In (l ± (6) 

where gi is the degeneracy factor and the upper sign is for 
bosons and the lower sign is for fermions. Here Xi is the 
particle fugacity, related to particle chemical potential 
fXi = Tin Aj. 

The yield average and fluctuation is then given by: 


{N,)gc 

(AA^^>gc 


dlnZ, 

dXi 

02 \nZi 



Anp^dp 

(27r)3 


rii (Ep), 


( 7 ) 



Airp^dp 

(27r)3 


ni{Ep) (1 =F ni{Ep)). (8) 


and 


ni{Ep) 


1 


(9) 


v{Q) = (AQ2) / (iv,h) (3) 

(where Ach = A+ -|- A_) proposed in the past as a probe 
of the QGP formation jj. First results for v{Q) are also 
available from RHIG experiments pm 

In the SHM, the charged particle multiplicity is given 
by summing all final state (stable) charged particle mul¬ 
tiplicities. These can be computed by adding the direct 
yield and all resonance decay feed-downs. The total yield 
of a stable particle a is 

(Na) total = {Afa)GC + '^Bj^a{Nj)GC (4) 

where j labels resonances. is the probability 

(branching ratio) for the decay products of j to include 
a. The charged particle multiplicity is given by the sum 
of all charged stable particles. 

The net charge fluctuation is given by 


We note that Xi enters the partition function in Eq. ©■ 
Hence, the validity of Eqs. o and 0 depends on weather 
Eq. 0 can be used as a generating function for the prob¬ 
ability distribution of states. It is important to underline 
this as in a dynamical system the value of Xi is not de¬ 
termined solely in terms of entropy maximization, but is 
subject to chemical conditions prevailing in the system, 
and here importantly, includes effects related to chemical 
non-equilibrium. Where Eq. ini represents a generating 
function but the system is not in chemical equilibrium, 
the fugacity Xi , is not anymore a Lagrange multiplier but 
a parameter characterizing the quantum number density. 

In a scenario where freeze-out occurs as a break-up 
of a chemically equilibrated hadron gas, the fugacity of 
the hadron i is given by the product of the fugacities of 
conserved quantum numbers. 

A^ = X^K-^xB , Af = (Ar)-\ (10) 


i 


where g, q is the number of light anti-quarks and quarks, 
respectively and s, s is the number of strange anti-quarks 
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and quarks, respectively and I 3 is the isospin. This for¬ 
mula implies that the fugacity for the antiparticle is, in 
full chemical equilibrium, the inverse of the fugacity for 
the particle, and the fugacity for a hadron carrying van¬ 
ishing conserved quantum numbers is 1 . 

In our approach, we do not assume that that the chem¬ 
ical equilibrium is reached nil. Hence Eq. no 
longer applies. The deviation from chemical equilibrium 
can be parametrized by a phase space occupancy factor 
7 q (for u, u, d, d in hadrons) and 7 ^ (for s and s). In this 
chemical nonequilihrium case the fugacity becomes 

A, = Ar7r^r^ (n) 

where A®'^ is given by Ea. ltTIIIl (Note that 7 , = 07 ). 

A system undergoing collective expansion is unlikely 
to be in chemical equilibrium, since collective expansion 
and cooling will make it impossible for endothermic and 
exothermic reactions, or for creation and destruction re¬ 
actions of a rare particle, to be balanced. However, since 
inelastic collisions have in general a slower relaxation 
time than elastic ones, an approximately perfect fluid can 
still have 7 7 ^ I (it’s evolution will be a non-trivial func¬ 
tion of time, since 7 does not commute with the Hamilto¬ 
nian). Furthermore, light quark chemical nonequilibrium 
is well motivated in a scenario where an entropy rich de- 
confined state quickly hadronizes |^. In this scenario, 
mismatch of entropies between the two phases requires 
Iq > 1 . 

Despite the lack of equilibrium and entropy maximiza¬ 
tion w.r.t. conserved quantum numbers, we will argue 
that the Eqs. |H1 and [T] apply in such a situation, with 
7 s contributing to the chemical potential via Eq. dJ. 
The validity of Eq. and 0 depend on the extent that 
Eq.© represents a probability generating function for 
the statistically hadronizing system. Within a statisti¬ 
cal hadronization scenario where hadrons are formed in 
proportion to their phase space weight given (not neces¬ 
sarily equilibrated) densities this is indeed the case 
provided the dynamics behind 7 does not generate addi¬ 
tional, non-statistical fluctuations. For an instance where 
the last issue is a concern, fluctuations of a quantum 
number produced mostly in initial-state processes (such 
as charm |^) will likely be dominated not by the 
statistical hadronization contribution but to fluctuations 
in initial abundance. 

Given that in the considered model non-equilibrium 
arises due to the ramd hadronization of the collectively 
expanding system [^, and since the observable charged 
particles are produced not in in the initial state but dur¬ 
ing the final break-up of a locally thermalized system, 
such non-statistical fluctuations should not be signihcant 
for the observable we are considering. Similarly, as we 
have argued in the previous section, initial-state volume 
fluctuations give a negligible contribution to the observ¬ 
able under consideration. 

However, it is possible that additional sources of irre¬ 
ducible two-particle correlations and fluctuations could 
arise near a phase transition. These effects go beyond 


the scope of this work. We will however argue that the 
applicability of our scenario, and the absence of further 
correlations can be tested by requiring that the same tem¬ 
perature and 7 s describe both the yields and the fluctu¬ 
ations of all soft hadronic observables. As we will show, 
this is a very stringent requirement. If it turns out that 
a single set of T, A®'^ and 7 ^ and 7 ^ is capable of de¬ 
scribing all yields and fluctuations, then it certainly is 
a strong indication that Eq. 0 can be interpreted as a 
generating function of the probabilities. The goal of this 
paper is then to hnd a way to experimentally determine 
the additional parameter 7 ^ which can be then used to 
compare the SHM calculation of yields and fluctuations 
to the experimental measurements. 


IV. FLUCTUATIONS IN CHEMICAL 
NON-EQUILIBRIUM 

Chemical nonequilibrium is of a particular interest 
since it can result in a large pion fugacity which influ¬ 
ences fluctuations much more severely than the yields. 

If 7 g becomes large enough so that A^ approaches 
gm„/T^ then the pion yield and the fluctuations behave 
like (c.f. Eas. (l7l8ll l 

lirn(7V) oc e“^, liin(AiV)^ oc e“^. (12) 

where e = 1 — \T^e~™"" 1"^. The fluctuation grows much 
faster than the yield as mentioned above. 

Some studies of yield ratios have indeed found the 
value of 7 q that can potentially make e small [s^, 

1^ . However, other studies of yield ratios concluded 
that 7 g is not necessarily large due to the fact parameters 
in such fits are highly correlated. In this case, adjusting 
other parameters such as the temperature can accom¬ 
modate current data without having 7 ^ 7 ^ 1 , but with 
much reduced statistical significance. Since such conflict 
is common when only the yields are considered, it be¬ 
comes necessary to study fluctuations as an additional 
constraint to determine the occupation factor jq more 
convincingly. 

We now discuss our specific analysis results. We used 
the public domain SHM suite of programs SHARE |^. 
expanded to include the fluctuations We evalu¬ 

ate yields and fluctuations, allowing for production of 
hadron resonances, their decay, and a possible absence 
of chemical equilibrium. In the rest of this paper, we 
set Af®' = I,A®'i = = 1.05 and A®'! = 1.027 in 

accordance with |30|. However, the two observables we 
consider, the net charge fluctuations and the A/K“ parti¬ 
cle yield ratio, are nearly independent of these quantities 
as will be shown below. 

Fig- d shows the variation in v{Q) as a function of jg 
for T = 140,170 MeV. The solid lines show v{Q) includ¬ 
ing the resonance decays, dot-dashed lines comprise only 
the direct effect of pion fluctuations. As the temperature 
increases (solid lines from top to bottom) the number of 
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FIG. 1: (Color online) v{Q) as function of 7 ^ (solid lines). 
Dot-dashed lines, no resonance decays; dashed lines, Boltz¬ 
mann fluctuations. Ellipses (blue) indicate the expected 
result areas for the equilibrium ( 7 ^ = 1 , solid) and non¬ 
equilibrium (75 yf 1 , dashed) models. 


resonances increases. This in turn increases the unlike- 
sign charge correlations and hence reverses the tempera¬ 
ture dependence of the pure pion case (dot-dashed lines). 
The short dashed lines show results for Boltzmann statis¬ 
tics. Boltzmann charge fluctuations are nearly constant 
as function of 7 ^ and primarily depend on chemical mix 
of the directly produced and secondary decay particles, 
which dominantly depend on the temperature T. The 
solid and dot-dashed lines in Fig. Q] terminate when the 
fluctuations start to diverge as in Eq. m- 

To determine both T and 7 ^ values we require an ad¬ 
ditional observable. In this work, we choose the yield 
ratio A/K~. This ratio depends linearly on 7 ^, and 
is nearly independent of and 7 s as A = (sdu) and 
K~ = (su). In Fig. El we show how the relative yield de¬ 
pends on 7 g and T. The A yield we wish to consider does 
not include weak decay feed from S but it includes the 
electromagnetic decay of and the strong decays. K~ 
excludes feed-down from (j), but includes K* and higher 
resonances. It is important to exclude the S and (j) cas¬ 
cading in order to eliminate the dependence on 7 s and 
A®'^. Fortunately, this is experimentally feasible. 

A similar ratio, which is experimentally easier to cor¬ 
rect for, is also dependent on temperature and 7 ^ 
only. See |37l| for the equivalent discussion in terms of 
E/cj). 

We now combine results in Figs.^andElinto our main 
result Fig. El Every point in this plane of v{Q) and A/K 
corresponds to a specific set of T and 7 ^ as indicated by 
the grid. Note that some domains in this plane are not 
allowed since they lie in the region where the (gener¬ 
ating, GC) partition function cannot be defined. The 
two highlighted regions indicate the expected chemical 
equilibrium (solid line ellipse at small v(Q), correspond¬ 



FIG. 2: (Golor online) Particle yield ratio A/A'“ as a func¬ 
tion of T (right panel) and 7 ^ (left panel) The A yield does 
not include H —> A and the K~ yield is without the contri¬ 
bution of ()) —> decays. Ellipses (blue) indicate the 

expected result areas for the equilibrium ( 7 ^ = 1 , solid) and 
non-equilibrium ( 7 ^ 7 ^ 1 , dashed) models. 

ing to 7 g = I and T = 170 MeV) and nonequilibrium 
parameter domains (dashed line ellipse at larger v{Q), 
corresponding to jq = 1.62 and T = 140 MeV). When 
particle yields and fluctuations are considered, the sepa¬ 
ration of these two domains confirms that we have found 
a sensitive method to determine both 7 ^ and T. The re¬ 
sults of having two extreme values, 7 s = 1 and 7 s = 2.5, 
are also shown in Fig. El The 7 s values corresponds to 
the equilibrium |35| and non-equilibrium [24| best fits. 
Their difference, as seen in Fig. El is small and well be¬ 
low the experimental error. The largest remaining sys¬ 
tematic deviation is due to the baryon chemical potential 
g/is/ST _ j^)g contribution to v{Q) is negligible, but 
this is not true for the case of A/K~. Generally the value 
of A®'^ is well determined by baryon to antibaryon yield 
ratios in a model independent way. 

To transform the diagram in Fig. El (or E/cj) in |^) 
to an equivalent result applicable to lower reaction en¬ 
ergy where A®'^ is greater, one has to allow for this 
change: We note that A/K~ oc (A®'^)^, and thus we need 
to multiply the axis in Figs. El and El by (A®‘i)^/1.05^. 
One can actually use the A/K ratio in this. Since 
AK^/AK~ oc (A®'^)®, the axis rescaling would be done 
with {AK'^/AK~y^‘^/l.Ob^ (A, AT corrected for 5 and (f) 
feed-down). 


V. ISSUES RELATED TO DETECTOR 
ACCEPTANCE 

The main phenomenological issue that prevents the 
straight-forward extraction of parameters from graphs 
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v(Q) 


FIG. 3: (Color online) Particle ratio KjK~ and particle fluc¬ 
tuation v(Q) plane: a point in plane corresponds to a set 
of values ')q,T. Black Lines correspond to results at fixed 
T = 200 (top), 170, 140, 100 MeV (bottom). The red dashed 
lines are for 7 , = 0.8, 1,1.4, 1.6, 1.8 from left to right. Thick 
lines correspond to 7 s = 2.5, thin lines correspond to 7 s = 1. 
Ellipses (blue) indicate the expected result areas for the equi¬ 
librium ( 7 q = 1 , solid) and non-equilibrium ( 7 ^ 7 ^ 1 , dashed) 
models. 


such as Fig. |21 are effects relating to the detector ac¬ 
ceptance. First of all, it has long been known that v{Q) 
is not a “robust” observable, but in general depends on 
the detector’s kinematic (rapidity and pr) cuts. This dif- 
hculty, however, can be lessened via mixed event back¬ 
ground subtraction. It can be shown 0 that observables 
corrected this way are in certain limits “robust” w.r.t. 
kinematic cuts and detector response. 

We have discussed how to generalize the methods de¬ 
scribed in this paper to robust observables elsewhere 
HE HE HE’ hence will not dwell on this topic, be¬ 
yond noting that, while diagrams such as Fig. 13 need 
to be re-thought since dynamical observables generally 
also depend on the (average) system volume, the sen¬ 
sitivities of the fluctuation and yield observables to the 
statistical model parameters follow the pattern described 
by this paper. Hence, generalizing the methods described 
by this paper to dynamical observables (whether via hts, 
as was done in [sEl or three-dimensional diagrams), is not 
a difficult task. 

An issue that needs to be addressed separately, how¬ 
ever, is the acceptance dependence of particle correla¬ 
tions. If the detector’s pseudo-rapidity coverage is too 
large, than the small volume assumption required for 
the Grand-Canonical ensemble becomes untenable, and 
long-range correlations (such as global conservation laws) 
can modify fluctuations. If the detector’s pseudo-rapidity 
coverage is too small, correlations due to resonance de¬ 
cays acquire a rapidity-dependent correction (which is 
not eliminated by mixed-event subtraction since it cor¬ 


rects two-particle correlations). We will address these 
issues in the next sub-sections. 


A. Influence of conservation laws on fluctuations 


If the detector can capture the full phase space of the 
system than, barring dramatic departure from standard 
model physics, the net charge of the event can not fluc¬ 
tuate. More generally, if the phase space size of the de¬ 
tected system becomes comparable to the total system 
size, observables will not anymore be given by the Grand- 
Canonical ensemble. 

If the system is a fluid (or in general not in global equi¬ 
librium) no ensemble is expected to provide a good de¬ 
scription of fluctuations beyond the small volume Grand 
Canonical limit, since the observable region of phase 
space will include many locally equilibrated volume el¬ 
ements exchanging energy and quantum numbers via hy¬ 
drodynamic flow. While yields could still be approxi¬ 
mated by some ensemble, the long range correlations and 
global non-equilibrium should break all simple scaling of 
fluctuations with yields. 

Hence, the configuration space coverage needed for a 
statistical description needs to be appropriately small for 
the corrections to the GC ensemble to be kept under 
control. 

To investigate these corrections quantitatively, con¬ 
sider the Taylor-expansion of the entropy of the “reser¬ 
voir” : 


^(fVtot - N) 


S{Ntot) - N 


-N^ - 

2 dm 


dS 

'm 


Ntc 


Ntot 

+ ... 


(13) 


where Wot is the total number of particles in the reser¬ 
voir and the small subsystem, and N is the number of 
particles in the subsystem. The first and second terms 
result in the usual Grand-Canonical ensemble result 
through the identification of the equilibrium chemical po¬ 
tential p = —T{dS/dN). 

The third term gives the first correction; The Grand- 
Canonical ensemble is therefore a valid approximation 
when 


Cgc 


im^^s/dN^ 

2 {dS/dN)N,^, 


(14) 


This quantity can be easily related to more common ther¬ 
modynamic quantities 


1 T (N) 

CGC = ^^^kvtot (15) 

where (N) is the average multiplicity of the observed vol¬ 
ume and kvtot is the susceptibility of the total volume. 
For the relativistic ideal gas, this is given by 


Cgc 


V 




(16) 
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and, as shown in section UH 

V At/ 

Vtot (At/) tot 

where Ar/ is the detector’s (pseudo)rapidity coverage and 
{Ar])tot is the system’s rapidity interval. 

Thus, we discover that the larger the susceptibility 
is, the smaller the system size V/Vtot has to be for the 
Grand-Canonical limit to hold. 

In fact, the physics determining the departure from 
this limit is precisely the same as the physics determining 
the divergence of fluctuations within an over-saturated 
pion gas. This is unsurprising, since over-saturation is 
argued for as a signature of a phase transition, and in 
finite systems undergoing phase transitions it is the finite 
size of the system that gives a cut-off for fluctuations. 

The pion chemical potential of the system created at 
RHIC, however, is kept below divergence, so it is hoped 
that one unit of rapidity, corresponding to V/( 2 Vtot) ^ 
7%, provides a safe limit for the Grand Canonical en¬ 
semble. In such a small rapidity interval, however, corre¬ 
lations due to resonances need to be suitably accounted 
for. The next sub-section shows how to do that. 


B. Disappearance of resonance correlations at 
small Ar/ 

If charge fluctuations are calculated after all resonances 
have decayed, then Eq. 0 becomes 

((AQ)2) = ((AiV+)2) -p ((AA_)2) - 2 {AN+AN_} 

(17) 

where the last term accounts for unlike-sign charge corre¬ 
lations coming from the decay of neutral resonances. For 
a conserved charge, and full acceptance of all resonances, 
this expression is equivalent to Eq. ®, with the correla¬ 
tion term exactly balancing out the amplification of res¬ 
onance abundance fluctuations through the greater mul¬ 
tiplicity of resonance decay products, within a hadron 
gas the correlation term will be given by decays of the 
resonance j into N+ and A_ 

{AN+AN_)=Y,b,^+-{Nj) (18) 

3 

while the fluctuation of each stable N± has to be ali¬ 
mented by contributions to it from resonance decays [l| 

((AiV±)2) = ^((AiV±)2)^+ (19) 

i 

+ 6,^,(l - 6,^,) (A,) + ((AA,)2) j 

For a finite acceptance window in general not all reso¬ 
nances produced can be reconstructed, even if the effi¬ 
ciency of the detector were 100%. Hence these contribu¬ 
tions must be weighted with acceptance weight factors. 


and this applies here in particular to the limited rapidity 
acceptance. For a neutral resonance j decaying into 
positive particles and n_ negative particles, three such 
coefficients are needed: 

Two will be the fractions of the positively charged and 
the negatively charged decay products which land in the 
acceptance window, and the third will give the fraction 
of the H— pairs that will land in the window. These coef¬ 
ficients will modify the branching ratios bj^i in Eq. ra 
and &j—_in Eq. m- 

If boost-invariance is a good symmetry, the first two co¬ 
efficients can be fixed to unity, since particles coming out 
of the acceptance region are exactly balanced by particles 
coming in. However, this is not true for the number of 
detectable pairs. If the resonance is out of the detector’s 
acceptance window it is impossible for all of it’s decay 
products to be in a window. Hence, Eq. 03 will have to 
include a term giving the percentage of resonances whose 
decay products are both within the detector’s acceptance 
region. 

((AQ)2) = ((AA+)2) + ((AA_)2) 

- 2RFiT,Ay){AN+AN_) (20) 


The dependence of the observed fluctuations on Rp is 
shown in Fig. 3] left panel. 

We note two effects not considered here and believed 
to be unimportant: 

1 ) the rescattering after formation is unlikely to alter 
Rf, since the typical momentum exchange in each col¬ 
lision the exchanged momentum (q) ~ T/3 tends to be 
considerably softer than what is required to bring parti¬ 
cles outside the acceptance region (in most decays, the 
characteristic momentum of the decay products in a reso¬ 
nance’s rest frame p* tends to be significantly larger than 
this value); 

2) The higher-momentum pseudo-elastic “regeneration” 
processes, where detectable resonances would be created, 
are also unlikely to modify Rf since, by local thermal 
equilibrium, two particles coming into the acceptance re¬ 
gion through kinematically allowed pseudo-elastic inter¬ 
actions will be balanced out by two particles originally in 
the acceptance region which come out as a result of the 
re-interaction. 

Thus, a measurement of fluctuations can still be relied 
upon to gauge the number of resonances present at chem¬ 
ical freeze-out. This underscores the importance of fluc¬ 
tuations as a probe for freeze-out dynamics. 

We now obtain Rp for a azimuthally symmetric perfect 
detector having a pseudo-rapidity coverage Ar/. We shall 
follow the formalism in to relate the resonance’s rest 
frame (denoted by *) to the lab frame. 

For both particles -I- and — to be within the detector’s 
acceptance region, —Ar//2 < 1 /+,!/- < Ar//2 where 


»?± = 2 log 


/Ej - - pl± 


tjEl -\-Pl±^ 


= In 


cot 



(21) 
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If all angular dependence in the resonance’s decay ma¬ 
trix elements is neglected (a valid approximation if 
many resonances are produced, with an approximately 
azimuthally invariant distribution)the fraction of de¬ 
tectable H— pairs will then be simply given by a phase 
space integral 


^+-{i1r,Ptr) 


(Pp* 


n 



( 22 ) 


where: 


An 


Ap 


An 


1 

1 

1 

+ 

_1 

0 

r ^ ~\ 

0 


0 


and the function 0 (z) is the usual step function 

0(z) = 0 z < 0 
0(z) = 1 z > 0 

Now, for two body decays this reduces to 

^+-{vr,Ptr) = ^ j ^ ®+- (23) 

while for three body decays we use the Monte-Carlo rou¬ 
tine MAMBO fo generate points in phase space. 

To calculate r]+ and ry_ from the resonance rest frame 
kinematic variables we Lorentz-transform to the lab 
frame, and get 


Thus, A 77 should be as small as possible, statistics per¬ 
mitting, due to the not easily controllable corrections 
described in section EH A subsequent SHM analy¬ 
sis of the experimental data can than calculate Rp for 
each resonance decay important for charge fluctuations. 
Hence, a v{Q) , properly corrected for experimental ac¬ 
ceptance, can be computed from SHM parameters via 
Eqs.lO and I120|) , and fed into Fig. Hand similar figures 
or fits [H> l37l l38l| . The computational tools needed to 
perform such an analysis have been published separately 
as open-source software |^ . 

It is important to underline that to perform this anal- 
ys4 i it is not necessary to understand the full freeze-out 
dynamics of the fireball (local temperature, flow field, 
hadronization hypersurface). It is enough to have a sen¬ 
sible parametrization of b~^ in terms of particle mass. 
This function is commonl y o btained from particle spec¬ 
tra at thermal freeze-out |^ . and is approximately lin¬ 
ear in particle mass. The question is weather we can 
extrapolate to chemical freeze-out conditions with 
enough precision in a model-independent way. The rela¬ 
tively mild dependence of Rp on b~^, together with the 
fact that hadronic re-interaction decreases the tempera¬ 
ture and increases the flow and the high viscosity of the 
hadron gas makes us confident that we can do it. 


VI. SUMMARY AND CONCLUSIONS 


PL± = ±pI± -k 
PT± = iPT± + 


Plr 

rriR 

Ptr 

ruR 



P*-PR A 
ER-^rruRj 



P*-PR \ 
Er -k mR J 


(24) 

(25) 


To get an over-all fraction of accepted resonances which 
will enter Eq. , one has to convolute Eq. 113) with a 
resonance distribution function in momentum space 


poo pAr){2 

Rf= dpTR / dr]RP{riR,pTR)^+-iriR,pTR) 

Jo J-Ari/2 

(26) 

where P{riR,pTR) is a suitable distribution function for 
resonances normalized to unity. A suitable function in 
the low energy region at mid-rapidity is 


P{vr,Ptr) 


^PR Im dmTRmfj^e-'>^'rR 


(27) 


We have performed this integral using a Monte-Carlo 
method. The result is shown in the right panel of Fig. 
01 We note that the most abundant resonance decays 
for charge fluctuations do not depend strongly on the 
inverse slope parameter b~^-. Going from 6 “^ = 200 MeV 
to 5“^ = 300 MeV while staying in the same rapidity bin 
changes the p —> tttt correction by at most 5 %, and the 
less abundant but more sensitive 77 —> 7 r''" 7 r“ 7 r*^ correction 
by no more than 20 %. 


We have studied in this work how a simultaneous mea¬ 
surement of charge fluctuations and a ratio such as K/K~ 
can differentiate between chemical equilibrium and non¬ 
equilibrium freeze-out, and to constrain the magnitude 
of the deviation from equilibrium as well as the freeze- 
out temperature. Our results show that it is possible to 
distinguish the chemical equilibrium freeze-out condition 
7 q = I [33 with r = 170 MeV [s^) from the chemical 
non-equilibrium freeze-out condition 7 ^ = 1.6 1 ^ . 

This is mainly due to the increase in the fluctuations in¬ 
herent to an oversaturated Bose gas, see Eq. m- 

We have further discussed the dependence of two- 
particle correlations on the detector acceptance region, 
and have shown that it can be calculated to a reason¬ 
able precision in a model-independent way. The “right” 
experimental detector acceptance for a detailed study of 
fluctuations, therefore, is one that is appropriately small 
yet sizable to ensure the appropriate ensemble under 
study is Grand-Canonical, provided that acceptance cor¬ 
rections to resonance decays are properly taken into ac¬ 
count using the methods described in section lV^ Quan¬ 
titative corrections to Grand Canonical yield/fluctuation 
relations for the best fit parameters can be estimated 
quantitatively via Eq. d 

Provided the detector acceptance region for a given 
fluctuation measurement is published, Eq. d can be 
used to calculate a correction coefficient Rp to the 
{NpN-) correlation for each decay of a neutral reso¬ 
nance. Using a calculated Rp for each resonance de- 

















FIG. 4: (Color online) Left: Sensitivity of the charge fluctuation measure on Rf, the fraction of resonance decay products which 
remains in the detector acceptance window (c./. Eq. 12()ll i. Thin black lines denote T = 170 MeV, thick red lines T = 140 MeV. 
Right: Acceptance fraction for different resonance decays as a function of the inverse slope b (c.f. Eq. iHZl) and the detector 
pseudo-rapidity acceptance Arj (NB: rj in this context means the pseudo-rapidity. Not to be confused with the decay of the rj 
particle, shown on the right panel of the image). Acceptance regions of A?; = 6,4, 2,1, 0.5, 0.1 are considered, top to bottom in 
descending order 


cay, together with the statistical model parameters, the 
charge fluctuation variable v{Q) can be calculated from 
Eqs.0 and (I20II . This v{Q) will still retain the sensitiv¬ 
ities to temperature and jq demonstrated in section HVl 
since 7 ^ impacts the primordial fluctuation terms rather 
than the correlation. It can therefore be used, together 
with a measurement such as K/K~ as in Fig. |21 or within 
a fit as in 113,13, to test the validity of the statistical 
model, unambiguously constrain its parameters, and dif¬ 
ferentiate between the high-temperature equilibrium and 
supercooled over-saturated freeze-out scenarios. 

It is our intent to perform a complete data analysis as 
outlined here, including consideration of acceptance cor¬ 
rections and of resonance decays, once final RHIC fluc¬ 
tuation data becomes available. 
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